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Abstract 

In the context of scalar-tensor models of dark energy and inflation, the dy- 
namics of vacuum scalar-tensor cosmology are analysed without specifying the 
coupling function or the scalar field potential. A conformal transformation to the 
Einstein frame is used and the dynamics of general relativity with a minimally cou- 
pled scalar field are derived for a generic potential. It is shown that the dynamics 
are non-chaotic, thus settling an existing debate. 
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1 Introduction 



Recent cosmological observations have established that the universe is very close to being 
spatially flat, corresponding to vanishing curvature index K in the Friedmann-Lemaitre- 
Robert son- Walker ( "FLRW" ) line element 



in comoving coordinates (t,r,9,tp). We use units in which the speed of light and New- 
ton's constant assume the value unity, the metric signature is —,+,+, +, and the other 
notations follow those of Ref. pQ. If K were exactly zero, it would not be possible to 
establish it observationally, due to unavoidable experimental errors. Notwithstanding 
this, the cosmic microwave background experiments [2] measuring a total energy density 
p of the universe close to the critical density p c = |^ (where H = a /a is the Hubble 
parameter and an overdot denotes differentiation with respect to the comoving time t) 
can be regarded as a reassuring validation of the theoretical prediction that the universe 
was taken extremely close to the K = state by inflation early in its history PJ. 

Another surprising discovery, obtained with the study of type la supernovae at high 
redshifts jlj, is that the cosmic expansion is accelerated. In the context of Einstein's 
gravity this fact is explained by postulating that the pressure P of the cosmic fluid 
satisfies P < — p/3. In fact, the best fit to the observational data requires even more 
exotic properties for the dark energy, the fluid that accounts for 70% of the energy 
density of the universe in Einstein's theory; if w = P/p denotes the effective equation 
of state parameter of the dark energy, there is marginal evidence that w < — 1 [5 . This 
range of values of the parameter w corresponds to a Hubble parameter that is increasing 
with time (superacceleration): H > according to the Friedmann equation 



where k = 8irG and G is Newton's constant. Most dark energy models are based 
on scalar fields and, if the universe really superaccelerates, models based on general 
relativity with a canonical, minimally coupled, scalar field are unviable. In fact, the 
energy density and pressure of such a scalar field are 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



+ r 2 (d6 2 + sin 2 9 dcp 2 ) 



(1.1) 



1 - Kr 2 



H = --(P + p) 



(1.2) 
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observational data requires a dynamical form of dark energy with an equation of state 
parameter changing with redshift, which would definitely exclude the cosmo logical con- 
stant as an explanation for dark energy (such a model is anyway unappealing because 
of the cosmo logical constant problem and the cosmic coincidence problem jHlIZI)- For 
this reason, alternative models in which the universe can superaccelerate have been con- 
sidered, including phantom fields with negative kinetic energy jS], scalar fields coupled 
non- minimally to the curvature [HI CHI; or alternative gravity theories. In this paper, we 
consider scalar-tensor extensions of general relativity. The theories in this class [TT) ITS] 
generalize Brans and Dicke's ^3] theory, exhibit many features in common with string/M 
theory ^3], are the arena for extended ^Hj and hyperextended [H)] inflationary scenar- 
ios, and are widely used in cosmology [T71 [TH] . Scalar-tensor gravity is used to model 
inflation in the early universe or dark energy in the present late time era of the universe. 
Both phantom cosmology and the theory of a non-minimally coupled scalar field, capable 
of producing the superacceleration phenomenon, can be seen as scalar-tensor theories. 

Because of the wide use of scalar-tensor gravity, motivated by the belief that the 
latter may be more fundamental than Einstein's gravity, one would like to have a clear 
and general picture of scalar-tensor cosmology. In general relativity the dynamics of a 
particular (dark energy or inflationary) scenario is usually determined once the poten- 
tial is fixed; in scalar-tensor gravity one also has free coupling functions adding extra 
degrees of freedom. However, it would be desirable to understand the dynamics with as 
much generality as possible without fixing these functions. This is what we set out to 
do in this paper and, on the basis of the observational data and of theoretical prejudice, 
we restrict ourselves to considering a spatially flat (K = 0) universe. In particular, we 
derive conclusions about the dynamics based on general assumptions about the poten- 
tial (e.g., monotonic, etc.). To the best of our knowledge, this classification has not 
been explored even in the relatively simple case of general relativity with a minimally 
coupled scalar field, in which the potential V(4>) is the only unknown function (e.g., in 
inflationary scenarios). Moreover, the geometry of the phase space is rarely discussed, 
even in scenarios based on general relativity with specific choices of the potential V(0), 
for which phase space studies exist in the literature (see Refs. JH]-[22] for Brans-Dicke 
theory) — usually only projections of the phase space are considered. The role of chaos 
in cosmology has received much attention since the early work on anisotropic universes 
approaching an initial singularity (mixmaster universes) |2*3j . 

A particularly interesting aspect of the cosmological dynamics is the presence or 
absence of chaos in a K = FLRW universe with a single scalar field; this has been the 
subject of debate. While the presence of chaos in the dynamics has been suggested on the 
basis of numerical studies and of a Painleve analysis [231122], it has been stated explicitly 
that chaos can not occur in the case of a single, minimally or non-minimally coupled, 
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scalar field • This has led to regarding inflation in the early universe as playing a new 
role: by taking the universe extremely close to a spatially flat one, primordial inflation 
eliminates the possibility of chaos from the dynamics. This statement is based on an 
analytic study of the phase space and its dimensional reduction to a two-dimensional 
surface [26, 27). Since this dimensional reduction can be generalized to any scalar- 
tensor K = cosmology with a single scalar field, it would appear that chaos is absent. 
However, the statement is based on the two-dimensional nature of the phase space 
but the known theorems (e.g., [2E]) do not apply to a multi-sheeted, non-compact and 
non-connected phase space such as the one of scalar-tensor cosmology. Therefore, the 
problem of the presence or absence of chaos remains open. We show in the following 
that, while determining the dynamics for general potentials, one is also able to exclude 
the presence of chaos in the dynamics of scalar-tensor K = cosmologies. 

In section 2 we approach the problem of determining the dynamics in the context of 
(spatially flat) scalar-tensor cosmology in the Jordan frame, while in section 3 we reduce 
it to the equivalent problem in the context of Einstein's theory by means of a conformal 
transformation to the Einstein frame. In section 4 we discuss the dynamics of the scale 
factor and of the scalar field for general potentials in general relativity, and then the 
conformal transformation is used to map back the results to the Jordan frame. Section 5 
contains a discussion and the conclusions. 



2 Scalar-tensor cosmology 



We consider the scalar-tensor class of theories jTTj described by the general action 

1 



S 



16n 



d X\/—g 



^g ab V a <pV b <p - V(<f>) 



(2.1) 



We do not add to the action terms describing "ordinary" matter (as opposed to scalar 
field 0) because we are interested in situations in which the cosmic dynamics are domi- 
nated by the scalar field, e.g., during early inflation or during the late time dark energy 
era. The field equations are 



R 



ab 



QabR 



u{<j>) 



V a 0V t 



1 



■ g ab V c 0V c ^ 



+ - (V a V fe - g a bO(, 



V{<j>) 



-g a b 



(2.2) 



2lu(<P) + 3 



0^-2K(0)-^V c 0V c 



(2.3) 
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where V a is the covariant derivative operator associated with the metric g a b, and □ 
g V V&. In a K = FLRW metric these equations assume the form 



H = —H — + 



, u{<t>) U\ v{<j>) 



6 



+ 



H 



1 



+ 13H + 



UJ 



2uj + 3 



<t>) 20 (2u(<f>) + 3) 
1 



<t>V -2V + J 



2uj + 3 



(2V - <f>V) 



(2.4) 

, (2.5) 
(2.6) 



where a prime denotes differentiation with respect to 0. Eqs. ()2.4|) - (j2.6j) are respectively, 
the Hamiltonian constraint H 2 = Kp^/3, H = —k (p^ + P^) /2, and the Klein-Gordon 
equation for 0, where p^ and P^ are the effective energy density and pressure of the 
scalar. Only two equations in the set (|2.4)) - (j2.6j) are independent. We choose H and 
as dynamical variables (other works in the literature choose different variables for 
ease of manipulation in special cases, but the physical meaning of their results is some- 
what obscured). To analyze the phase space it is convenient to regard the Hamiltonian 
constraint (|2.4J) as the algebraic equation for 



u;0 2 



with roots 



(H, ' 



+ (<f)V - 6H 2 <j) 2 ) = 



3H<j) ± ^3H 2 <P 2 (2uj + 3) - uj<pV 



(2.7) 



(2-8) 



Unless uj = 0, the phase space is a two-dimensional surface S = S + U X~ composed of 
two sheets E corresponding to the positive or negative sign in eq. (|2.8|) . embedded in 
the three-dimensional space (if, 0, 0). In actual fact, only the region > is of physical 
interest because it corresponds to positive gravitational coupling, and there are physical 
constraints on the variability of the effective gravitational coupling 



[UJ 



2) 



(2uj + 3) 



(2.9) 



(see, e.g., Refs. H2HB|). 

In general, there is a dynamically forbidden region 



F={{H,<j>)\ 3H 2 (f) 2 [2u(<f>) + 3] - cj(0)0K(0) < 0} 



(2-10) 
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corresponding to a negative argument of the square root in eq. (|2.8|) . This forbidden 
region may be absent for specific choices of uj{4>) and V(<p), but is present in most 
scenarios considered in the literature (see Ref. [21] for further details). The boundary 
of the forbidden region is composed of the only points where the two sheets touch each 
other and <fi is single valued, i.e., the set 

B = \[H,<j>): 0=^| = S + nS-, (2.11) 



or 

r2 i2 



(H,(j)): 3H 2 (j) 2 (2uj + 3) =(f)ujV (2.12) 

on B, which defines implicitly a curve H(<j)) in the space (H, (p, </>). 

Having chosen H and <fi as dynamical variables, the equilibrium points of the system 
are necessarily de Sitter spaces with constant scalar field (H ,(j) ). If they exist, these 

can lie anywhere in the phase space Sfl |</> = o|. According to the dynamical equations 

fl2.4}) - (j2.fij) . there are two necessary and sufficient conditions for the existence of such de 
Sitter fixed points. 



r2 V(<j> ) 



600 

00^-2^ = 0. (2.14) 

This two-dimensional structure of the phase space was reported earlier in the literature 
for minimally [3U] and non-minimally coupled [HTJ E2j scalar fields. In Ref. it was 
argued that the reduction of the phase space to two dimensions implies that chaos is 
impossible, and therefore that primordial inflation taking the universe extremely close 
to an exactly spatially flat FLRW space in addition to solving the problems of standard 
Big Bang cosmology [Sj, has the effect of inhibiting chaos. However, this conclusion is 
based on the Poincare-Bendixson theory which applies to a two-dimensional phase space 
that is flat instead of curved, and to regions that are compact and connected. This is 
certainly not the case of the £ phase space which, in addiction to being double-sheeted 
and curved, extends to infinity in all directions and in general is not connected due to 
the presence of the forbidden region, which can consist of two or more separate "holes" 
in £ (see Ref. j^Zj for examples). 

At best, the argument of Ref. [2H] arguing against the presence of chaos for the special 
case of nonminimally coupled fields can be applied (and generalized to arbitrary scalar- 
tensor cosmologies described by eqs. ()2.4|) - ()2.6^ to compact regions of the phase space £ 
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as follows. Consider a compact region C in one of the sheets S + or £~, lying away from 
the boundary B of the forbidden region. Let (u, v) be smooth local coordinates covering 
C (if u and v do not cover the entire region C one can consider an atlas composed of two 
or more charts); 

u = u{H,4>) , (2.15) 



v = v{H,4>) , (2.16) 

is a smooth map from a region of IR 2 to the sheet X 1 * 1 , which is locally flat. Then, using 
the comoving time t as a parameter, 

,/ 



du ■ du ■ 
u = -—H+—(t> = 

dH dcj)^ 

(f)V'-2V du 
+ 20(2w + 3) dH 

dv ■ dv ■ 

v = m H + T4> <p = 

<f)V'-2V dv 
+ 2<j)(2uj + 3) dH 



+ 



UJ 



du 



(2H_du_ du\ 



2cj) 2 2(f) (2u + 3) J r dH \ <p dH d<p J 



+ 



LO 



2 20 (2oo + 3) 



, dv 

1m 



+ 



2H dv dv 



dH 



(2.17) 



(2.18) 



By using the expression (J2.8|) of <f> and the fact that u, v, and their derivatives only 
depend on H and 0, one can write eqs. (|2.17|) and (J2.18|) in the form of the autonomous 
system of first order ordinary differential equations for (u(t),v(t)) 



u 



f{u,v) 



(2.19) 



g(u,v) 



(2.20) 



where 



/ (u, v) 



LO 



+ 



LO 



(2lo + 3) 



3^0 ± v/3iJ 2 2 (2u + 3) - uu<pV~ 



LO 



du 
dH 



+ 



+ 



3#0 ± ^3H 2 <f) 2 {2u + 3) - u<f)V (2H du du 



LO 



\ dH 



+ TTT 



W -IV du 
20 {2u + 3) dH 



(2.21) 
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and 



g(u,v) 



OJ 



+ 



OJ 



20 2 20 (2u + 3) 



3F0 ± ^3H 2 <p 2 ( 2u j + 3) - u<f)V 



OJ 



dv 
dH 



+ 



3H(j) ± ^3H 2 (f) 2 {2u + 3) - u(f)V ( 2H dv dv 



OJ 



dH 



<f)V'-2V dv 



(2.22) 

+ 20(2^ + 3)^' (2 ' 23) 

The reduction of the dynamical system to a first order autonomous system with phase 
space consisting of a plane allows one to apply the standard Poincare-Bendixson theory 
which guarantees that compact and connected regions of E 1 * 1 (corresponding to compact 
and connected regions of R 2 ) are free of chaos. However, this argument only proves the 
statement of Ref. [251 f° r sucn compact regions but not for the entire phase space E. 
This larger phase space can be studied in the context of a more comprehensive analysis 
of the dynamics that we propose in the next two sections. 



3 Conformal mapping of the phase space 



It is well known that the Jordan frame action ()2.1|) of scalar-tensor gravity can be 
mapped to the Einstein Hilbert action by means of the conformal transformation 



9ab — ► 

and the scalar field redefinition 



9ab — Q 2 g a b, ^ 



given by 



(3.1) 



2u(<f>) + 3 
16V 



(3.2) 



(see Refs. |32[ 13*3*1 El for reviews). In terms of conformally rescaled quantities in the 
Einstein frame, which are denoted by a tilde, the action (|2.1j) assumes the form 



S 



d Xa 



R 



1 



16ttG 



r b Va0V fe <; 



U 



where 



V 



(3.3) 



(3.4) 
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The sign of the potential is preserved by the conformal transformation. By studying the 
dynamics in the Einstein frame and mapping back to the Jordan frame, one can study 
the dynamics and, in particular, the absence or presence of chaos. 

Before we proceed, let us study the conformal cousins in the Einstein frame of the 
de Sitter fixed points living in the Jordan frame. It is easily concluded that de Sitter 
equilibrium points in the Jordan frame are mapped into de Sitter equilibrium points in 
the Einstein frame, and vice-versa. In fact, by using eq. (|3.2j) and the usual range of 
values of the coupling function 2a; (0) + 3 > given by Solar System experiments and 
guaranteeing a positive-definite kinetic energy term for the scalar, one has d<p/dt = if 
and only if dcp/dt = 0. In addition, eqs. and ()3.1|) yield 

d~s 2 = tfds 2 = n 2 (-dt 2 + a 2 dx 2 ) = -dt 2 + ~a 2 dx 2 , (3.5) 

where dt = \f<pdt and a = \f$a. Hence, 



dcj) /2aH-3 3/2 



dt V 167T dt 



(3.6) 



Because we consider positive gravitational coupling > 0, the sign of difi/dt is the same 
as that of 0, and d<p/dt = if and only if dip/dt = 0, while 

dH _ 1 / 30 2 ■ H<j>\ 
dt 1^20 403/ 2 + 2^ 1 ' 

Therefore, a de Sitter fixed point with constant scalar field [H, 0j = (0, 0) in the Jordan 

frame corresponds to a de Sitter fixed point with constant scalar field (^dH/dt, dcp/dt^j = 

(0, 0) in the Einstein frame (this fact is related to the scale invariance property of the 
exponential function). It is not true, as is instead stated in Ref. [33] that the conformal 
transformation becomes singular at the equilibrium points. 

The condition (j2.14J) for the existence of the Jordan frame fixed point (Hq, 0o) trans- 
lates into the corresponding condition dU/d(f> = 0. In fact, 



dU dU dd> 16vr 



and eq. (|2~T4]) implies ^ 



2^(0) + 30 

0, whereas = 0o = constant is equivalent to 



(3.8) 



<t>0 

t>o = constant. Moreover, stability of the de Sitter fixed point in the Jordan frame 
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corresponds to stability of the corresponding de Sitter point in the Einstein frame, as 
shown in the following. Since 



~ _ 1 da (p H 

lYt = W 2 + 7^' (3 ' 9) 

perturbations 5H and 6(f) in the Jordan frame correspond to the Einstein frame pertur- 
bations 

6H = 5^3/2 + -JT ~ o j. 3/2 °$ ' 3 " 10 



H = JMM±1^, (3 .n) 

V 1D7T o 

to first order. If the Jordan frame fixed point (H , O ) is stable, perturbations 5H and 
5<p do not grow (or decay rapidly) and the same can be concluded for the Einstein frame 
perturbations 5H, S(f>; and vice-versa. 



4 Dynamics in the conformally rescaled world 



We finally proceed to study the dynamics in general relativity in the conformally rescaled 
world. For economy of notations we drop the tilde and in this section an overdot de- 
notes differentiation with respect to the Einstein frame comoving time. The Einstein- 
Friedmann equations are 



H 2 



H 



K 

6 



2 + 2C/(0) 



K 

2 



+ 3#0 + 



dU 



. 



(4.1) 



(4.2) 



(4.3) 



Since general relativity is a trivial case of scalar-tensor gravity, the structure of the 
phase space discussed in section 2 is still present, albeit some simplifications occur. The 
Hamiltonian constraint (|4.1|) can be seen as an algebraic equation for <p with roots 



:(H, 




3H 2 



K 



(4.4) 
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which makes evident the double-sheeted structure of the two-dimensional phase space 
S = S + U S~ embedded in the three-dimensional space (H, 0,0). In general there is a 
forbidden region 

T = {{HA): —<U{if,)} (4.5) 



with boundary 

B = s + nr = {(#,0,0) : — = £7(0), = 0}. (4.6) 

The de Sitter equilibrium points (i£ , 0o) are the only de Sitter spaces, while in general 
scalar-tensor cosmology de Sitter spaces with non-constant scalar field may be solutions, 
although they are not fixed points (cf. eq. (|2.5|1 ). Moreover, contrary to the case of 
section 2, all the de Sitter fixed points are forced to lie on the boundary B = E + n E~ 
between the upper and lower sheet, corresponding to = 0. Eqs. (|4.ip - ()4.3|) yield the 
two conditions for the existence of de Sitter spaces (H , O ) (there are only two conditions 
because only two of the eqs. (|4.1j) - (|4.3|) are independent), 

Hi = I U , (4.7) 

U' = , (4.8) 

where U = U (0 O ) and U' Q = ^\<f> - de Sitter fixed points exist if the potential £7(0) has 
maxima, minima, or inflection points at 1; 2 , . . . , n , . . . with £7'(0i) = £7'(0 2 ) = . . . = 
£7'(0 n ) = . . . = 0. One expects that if £7 has a maximum (respectively, a minimum) 

in (pi, the corresponding fixed point (y/ kU (0j) /3, 0jj will be unstable (respectively, 

stable) as we prove later. 

Eq. (|4.2j) allows one to conclude immediately that there are no limit cycles (periodic 
orbits) because H is always decreasing (apart from the fixed point solutions) and can 
not come back to a previous value during the evolution of the system. Each point of 
B where U' 7^ can be crossed by a trajectory only once, and in a definite direction 
specified by the sign of U' at that point. 

In the upper sheet with boundary removed £ + \£>, it is > and H < 0, while in 
Y,~\B it is < and H < 0, hence can not change sign away from the boundary 
B. On B = S + fl X - it is H — 0, = 0, = — £7(0), as follows from the dynamical 
equations. The tangent to the orbits of the solutions parametrized by the time t in the 
space (^H, 0, 0j is the vector 

f(t) = (iy(t),0(t),0(o) (4.9) 
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and, on £>, 



f(t)\ B = {0,0, -U') 



(4.10) 



On the boundary B, the tangent T can only be vertical and no motion along the H or 
directions can occur. This means that there can be no motion along the curve B and 
oortions of B can not be parts of orbits of solutions. At points of B where U' = 0, it is 

T =0 and there is no motion: these points are de Sitter fixed points. At points of 

B where U' > 0, the tangent T points downward along the negative 0-axis and an orbit 
can only go from the upper sheet S + to the lower sheet £~ by crossing the boundary 
B. If instead U' < at points of B, the tangent T to the orbit is pointing upward 
in the positive ^-direction and the orbit crosses from S~ to S + . This excludes the 
possibility that an orbit "bounces" on the boundary B back to the sheet it came from. 
The possibility is not excluded that, in certain potentials, once the orbit has crossed, it 
changes component of the vertical velocity <p and comes back to the boundary to change 
sheet again. From eqs. ()4.1|) - ([4.3|) it follows that 



H + 3H 2 = kU 



(4.11) 



Since H always decreases monotonically (except, of course, at the fixed points) either 
H tends to a finite limit (horizontal asymptote) Hq as t — > +oo, or else H — > — oo. (For 
a given potential U(4>), there may be orbits going to a finite Hq and other orbits going 
to — oo, depending on the initial conditions). 

If H — > H as t — > +oo, then H — > and eq. (|4.2j) implies that <fi — > 0, which implies 
that also the scalar field </> approaches a horizontal asymptote <f)Q. In this situation, 
eq. (j4.11|) yields 3Hq = kU (0o), i-e., we have a de Sitter fixed point satisfying eq. ()4.7|) . 
which can only exist if U'(<po) = 0. If U' ^ for all values of (i.e., if U((fi) is strictly 
monotonic), there is the possibility that H — > — oo instead, which corresponds to a Big 
Crunch in a finite time, but there is also the possibility that Hit) — > H asymptotically 
while <p(t) — > ±oo and both H, <fi — >• as t — ► +oo. An example is the exponential 
potential 



U (0) = A exp 
which gives power-law inflation [Hj 



± 



167T 
P 



(P>1) 



a(t) 



a t p 



(4.12) 



(4.13) 




(4.14) 
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for which if = p/t — > while — > +00 and — > 0. 

We can now consider the situation in which if — > — 00 and consider times sufficiently 
large so that H < 0. To proceed we need to make some assumptions on the potential: we 
first assume that U is monotonic, so U' has definite sign. The first possibility is U' < 
for all 0; then the Klein-Gordon equation = — 3ff0 — U' implies that in the upper 
sheet X + , where > 0, it is > and therefore — ► +00; the point (if, 0) representing 
the universe in the phase space goes to infinity without oscillating. In the lower sheet 
E~, where < 0, either has a horizontal asymptote with — > 0o, or else — > —00. 
The first possibility is excluded because it would imply that — > and, according to 
eq. ()4.2jh also if — > 0, which contradicts the assumption that if — > —00; hence also in 
this case (if, 0) goes to infinity without oscillations. 

In the case U' > for all values of 0, take points in the upper sheet E + with > 0: 
then either 0(t) — > 0o (horizontal asymptote) or else 0(t) — > +00. If has a horizontal 
asymptote 0o, i.e., if — ► 0o, then 0(t) — ► and also H(t) — > 0: but this contradicts the 
assumption that H{t) — > —00. Then it must be 0(t) — > +cxd and the point (if, 0) goes 
to infinity without oscillations. If instead we pick a point in the lower sheet S~ with 
< 0, then = — 3ff0 — U' < and 0(t) goes to minus infinity without oscillations. 

We can then consider the trivial case in which U' = for all values of 0: this represent 
a cosmological constant and gives unviable early universe models in which inflation does 
not stop or dark energy models that are disfavoured (a scalar field is introduced to get 
away from the cosmological constant and its problems |7j); therefore, this is a case of 
mathematical interest that we include for completeness. In this case = — 3ff0, which 
integrates to = C/a 3 where C is a non-zero constant, and — > 00: the point (if , 0) 
tends to infinity without oscillations. 

Generally, we can consider a non-monotonic potential with a minimum, or possibly 
several maxima and minima. In this case the sign of U(4>) changes as evolves; this 
situation is best studied with a Ljapunov function. 

First, assume that U(4>) has a single absolute minimum Uq = £/(0o) attained at the 
single value 0o of the scalar field (examples are the widely used potentials U = m 2 2 /2 

or U — A0 4 ), then (ff , 0o) = \ kUo/3, 0oj is a fixed point. The function 

L(ff,0) = ^ + f/(0)-f/ o (4.15) 

is a Ljapunov function. In fact, £7(0) > V0 7^ O and L (if, 0) > V (if, 0) 7^ (ff , 0o). 
Moreover, L(H , O ) = since at the fixed point = 0. Along the orbits of the solutions 
we obtain, upon use of the Klein-Gordon equation, 

JL = (0 + U') = -3ff0 2 < (4.16) 
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when H > 0, i.e., for all expanding universes. This guarantees that (H ,<j) ) is an 
attractor with an attraction basin at least as wide as the H > region of S. 

The Ljapunov function L(H, 0) also allows us to immediately conclude that con- 
tracting de Sitter spaces with H = —a/ kUo/3 are unstable because H < and L > 
in a neighbourhood of this equilibrium point (it is well known that contracting de Sitter 
spaces are unstable — see, e.g., Ref. [33]). 

If U(4>) has an absolute minimum that is assumed at two (or more) different values 
of 0, say 0i,02, ... ,0n, ... then U(4>) > U((f>i) = U(4> 2 ) = ... = U(<p n ) = ... = 
Umin V0 {0i, 2 , ... , n , . . . } and U(4>) must assume maxima between the minima. 
In this case the phase space £ will contain the attraction basins of the stable fixed 
points corresponding to the minima of U(<fi) (the maxima corresponding to unstable fixed 
points), with a separatrix going through a fixed point between two adjacent attraction 
basins 

5 Discussion and conclusions 

Based on the reported marginal observational evidence for present superacceleration of 
the universe, which can not be explained by general relativity with a canonical scalar 
field, we consider spatially flat homogeneous and isotropic cosmologies in scalar-tensor 
gravity. The dynamics are conveniently studied by a conformal mapping of the Jordan 
frame scalar-tensor theory into the Einstein frame, in which gravity reduces to general 
relativity and the rescaled scalar field is minimally coupled to gravity (a non-minimal 
coupling of the scalar to ordinary matter is irrelevant here because we consider the scalar 
to be the dominant form of "matter" in the universe, and other material sources can be 
neglected). Apart from the original motivation to study scalar-tensor gravity by using 
the mathematical trick of the conformal transformation to the Einstein frame, the study 
of the phase space and of the dynamics of general relativity with a minimally coupled 
scalar field is interesting in itself. We emphasize that this work is limited to spatially 
flat (K = 0) FLRW models because they describe our universe according to the recent 
observations of the cosmic microwave background. However, from the dynamical point of 
view, the K = +1 cases are even more interesting j3Hj, and even spatially flat universes 
containing Yank-Mills fields exhibit chaotic oscillations of these fields |40j. 

The results obtained about the Einstein frame dynamics can be mapped back to the 
Jordan frame. The general picture obtained in the Einstein frame is that of a phase 
space that is free of chaos, with the orbits of the solutions converging to attractor points 
or going to infinity. Depending on the form of the scalar field potential, there can 
also be power-law attractor solutions, which are well known to exist in several scalar- 
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tensor gravity theories |17j . When mapped back to the Jordan frame, the non-chaotic 
Einstein frame dynamics translates into non-chaotic Jordan frame dynamics; when they 
exist, stable (unstable) equilibrium points are mapped into stable (unstable) equilibrium 
points, their attraction basins are deformed but they are still present, and the boundaries 
between different attractor points are well defined separatrices also in the Jordan frame. 
The possibility of such boundaries having fractal dimension, which would be a clear 
signature of chaos |28., is ruled out (note that fractal basin boundaries have instead 
been found in K ^ FLRW, or in anisotropic universes 36 ). 

The dynamics of spatially flat scalar-tensor cosmologies is thus well defined and 
chaos-free, due to the dimensional reduction of the phase space £ to two dimensions, 
as conjectured in Refs. [2U I2Z|- However, the proof of this conjecture is non-trivial 
due to the complicated structure of the phase space £ = S + US". Only spatially flat 
cosmologies enjoy the reduction of the phase space to two dimensions: the scale factor 
a(t) of the FLRW metric (jl.lj) only appears in the combination H = a/a and its first 
derivative H in the field equations with K = 0. If K ^ instead, terms of the form K/a 2 
will appear, spoiling the reduction of the phase space to two dimensions. The surface £ in 

the three-dimensional space (^H, 0, <pj separates the orbits of the solutions corresponding 

to K — +1 (located above £ + or below £~) from orbits corresponding to K — — 1 
(located between £ + and £~). This property was first realized in Ref. [37] for the special 
case of chaotic inflation with a massive scalar field in general relativity, and it corresponds 
to the impossibility of having dynamical transitions between different topologies of the 
spatial sections of a FLRW universe. Orbits of the solutions corresponding to K = ±1 
are free to move in the three dimensions (H, 0,0), where there is "enough room" for 
them to wind around each other, so chaos can occur, and it has indeed been reported in 
the literature 

Because inflation takes the early universe extremely close to, but not exactly to, a 
K = FLRW model, the orbits of the solutions can in principle depart slightly from the 
surface £, but they are extremely close to it. Then, the possibility of chaotic dynamics is 
not ruled out exactly, but is postponed for an extremely long time. Of course, if a second 
scalar field is present and significantly contributes to the dynamics of the universe, then 
the dimension of the phase space jumps up by two and chaos becomes possible, and is 
also reported in the literature (e.g., Ref. jSHl) — chaos in the dynamics of two mutually 
coupled scalar fields is an important element of reheating after inflation. 

To summarize, in addition to taking a standpoint on the possibility of chaos and 
generalizing the context of this debate to any K = scalar-tensor cosmology described 
by the action (|2.1|) . we have provided a general picture of the phase space and of the 
dynamics (see also Refs. [IH]-[22] for the special case of Brans-Dicke theory, the prototype 
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of scalar-tensor gravity). Special choices of the arbitrary functions iv((f>) and V((f>) will 
be the subject of future research. 
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